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1 INTRODUCTION 

It is well known that both stars and cold interstellar gas play 
an important role in the gravitational instability of galac- 
tic discs (e.g., Lin & Shu 1966; Jog and Solomon 1984a, b; 
Berlin & Romeo 1988, and references therein). The local 
stability criterion for two-component discs of stars and gas 
■ can be expressed in the same form as Toomre's (1964) sta- 
5^ bility criterion, Q > 1, provided that Q is redefined ap- 
, propriately. Berlin & Romeo (1988), Elmegreen (1995), Jog 
(1996), Rafikov (2001) and Shen & Lou (2003) calculated the 
effective Q parameter as a function of the radial velocity dis- 
persions and surface mass densities of the two components. 
Shu (1968), Romeo (1990, 1992, 1994) and Wiegert (2010) 
evaluated the stabilizing effect of disc thickness, which is 
usually neglected but significant. 

As galactic discs contain both stars and gas, the effec- 
tive Q parameter is clearly more accurate and useful than 
Toomre's Q. Berlin et al. (1989a, b) and Lowe et al. (1994) 
showed that the radial profile of this parameter has a large 
impact on the dynamics and evolution of spiral structure 
in galaxies. The results of those comprehensive analyses are 
discussed further in the books by Berlin & Lin (1996) and 
Berlin (2000). The effective Q parameter is also a useful di- 
agnostic for exploring the link between disc instability and 
star formation in galax;ies (e.g.. Hunter et al. 1998; Li et 
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ABSTRACT 

The Wang-Silk approximation, « Q^ars + Qgas : frequently used for estimating 
the effective Q parameter in two-component discs of stars and gas. Here we analyse 
tliis approximation in detail, and sfiow how its accuracy depends on the radial velocity 
dispersions and Toomre parameters of the two components. We then propose a much 
more accurate but still simple approximation for the effective Q parameter, which 
further takes into account the stabilizing effect of disc thickness. Our effective Q 
parameter is a natural generalization of Toomre's Q, and as such can be used in a 
wide variety of contexts, e.g. for predicting star formation thresholds in galaxies or for 
measuring the stability level of galactic discs at low and high redshifts. 
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al. 2005, 2006; Yang et al. 2007; Leroy et al. 2008). More 
applications and references are given below. 

Wang & Silk (1994) proposed a remarkably simple 
recipe for computing the effective Q parameter in the case of 
infinitesimally thin discs: Q^^ +Qg^, where and 

Qg are the stellar and gaseous Toomre parameters. Berlin 
(private communication) points out that such an approxi- 
mation was already used by him, before the 1990s, for il- 
lustrating the efficiency of a small amount of cold gas to 
destabilize a disc (see also Berlin 1996; Bertin & Lin 1996). 
Jog (1996) pointed out that the Wang-Silk approximation is 
invalid since it results from an incorrect analysis. In spite of 
that, the Wang-Silk approximation has been used in several 
important contexts: star formation (e.g., Martin & Kenni- 
cutt 2001; Boissier et al. 2003; Corbelh 2003; Wong 2009), 
galaxy formation and evolution (e.g., Immeli et al. 2004; 
Naab & Ostriker 2006; Kampakoglou & Silk 2007; Stringer 
& Benson 2007; Wetzstein et al. 2007; Foyle et al. 2008; 
Benson 2010), gravitational instability of clumpy discs at 
low and high redshifts (e.g., Bournaud & Elmegreen 2009; 
Burkert et al. 2010; Puech 2010), and others (e.g., Hitschfeld 
et al. 2009; Wong et al. 2009). 

In spite of such a burst of applications, there has been 
no attempt to assess how good the Wang-Silk approximation 
is. In this paper, we evaluate its accuracy by performing a 
rigorous comparative analysis (see Sect. 2.1). Besides, we in- 
troduce a new approximation for the effective Q parameter: 
simple, accurate and applicable to realistically thick discs 
(see Sects 2.2 and 2.3). We also show how to use our effec- 
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Figure 1. The parameter plane populated by nearby star- forming 
spirals. The galaxy data are from Leroy et al. (2008), Q* and 
Qg are the stellar and gaseous Toomre parameters, ct^ and erg 
are the radial velocity dispersions of the two components. The 
shaded part of the (s, q) plane is the two-phase region discussed 
in the text. The dispersion relation tj^(fc) has two minima inside 
this region, and one minimum outside it. The transition between 
the gaseous and stellar stability phases occurs for q = 1. This 
line intersects the boundaries of the two-phase region at {s,q) ~ 
(0.17, 1), where the stability threshold is Qgj^ ~ 1.4. 

live Q parameter for measuring the stability level of galactic 
discs, and why such a diagnostic is more predictive than the 
classical Toomre parameter (see Sect. 2.4). The conclusions 
of our paper are drawn in Sect. 3. 

2 APPROXIMATING THE EFFECTIVE Q 

2.1 Disc instability and 

the Wang-Silk approximation 

Wang & Silk (1994) investigated the link between star for- 
mation and disc instability in galaxies. They reconsidered 
the two- fluid dispersion relation of Jog & Solomon (1984a), 
which is valid for infinitesimally thin discs of stars and gas, 
and found that the effective Q parameter can be approxi- 
mated as follows: 



Qws Q* Qg 

where = kct^/ttGE* and Qg — nag/nGY^g are the stellar 
and gaseous Toomre parameters. This approximation is ap- 
pealing because it is as simple as the formula for the total 
resistance in a parallel circuit. To evaluate the accuracy of 
Eq. (1), we rewrite it as Qws = Q*/ QwS' where 

Qws = 1 + §^ • (2) 

The local stability criterion, Qws > 1, translates into Qj, > 
Qws- This is of the same form as the local stability criterion 



found by Berlin & Romeo (1988); see also Romeo (1985). 
Bertin & Romeo (1988) determined the stability threshold 
Qgjj^ numerically, starting from the same dispersion rela- 
tion as Wang & Silk (1994) but without introducing further 
approximations. In contrast to Qwsi Qbr depends on two 
parameters: o"g/cr* and Eg/E*. Since Eg/E* = (TgQ*/(7*Qg, 
we can easily express Qbr terms of 




We can then compare Q^c^[q) with Qbr('*i9) ^^'^ evaluate 
the accuracy of the Wang-Silk approximation as a function 
of s and q. 

Let us first see how spiral galaxies populate the (s, q) 
plane. We use the 12 nearby star-forming spirals analysed 
by Leroy et al. (2008), namely NGC 628, 2841, 3184, 3198, 
3351, 3521, 3627, 4736, 5055, 5194, 6946 and 7331. These are 
galaxies with sensitive and spatially resolved measurements 
of kinematics, gas surface density and stellar surface density 
across the entire optical disc. For each galaxy of this sam- 
ple, we compute the radial profiles s — s{R) and q — q{R), 
and hence the track left by the galaxy in the (s, q) plane. 
The result for the whole sample is shown in Fig. 1. The data 
span a range of two orders of magnitude in s and q, so we 
show them using a Log-Log plot. The typical value of s can 
be robustly estimated by computing the median of the data 
points along s, which is Smod — 0.27. This value is compa- 
rable to that found in the solar neighbourhood (s ~ 0.2; see 
Binney & Tremaine 2008, p. 497), but is much smaller than 
that expected in high-redshift star-forming galaxies (s ~ 1; 
e.g., Burkert et al. 2010; Krumholz & Burkert 2010). The 
median value of q, gmcd — 1.5, is close to unity. This sug- 
gests that, on average, stars and gas contribute equally to 
the gravitational instability of the disc. Similar values of q 
are found in the solar neighbourhood {q « 0.6; see Binney 
& Tremaine 2008, p. 497), and are also expected at high z 
{q ~ 1; e.g., Burkert et al. 2010; Krumholz & Burkert 2010; 
Tacconi et al. 2010) Q Last but not least, note that 20% of 
the data fall within the shaded part of the (s, q) plane. This 
is the 'two-phase region' of Bertin & Romeo (1988), here 
shown using our parametrization and logarithmic scaling. In 
this region, the contributions of stars and gas to the grav- 
itational instability of the disc peak at two different wave- 
lengths. If g < 1, then the gaseous peak is higher than the 
stellar one and gas will dominate the onset of gravitational 
instability. Vice versa, if g > 1, then stars will dominate. 
These are the gaseous and stellar stability 'phases' shown 
in Fig. 1. In the rest of the parameter plane, the dynamical 
responses of the two components are strongly coupled and 
peak at a single wavelength. More information is given in 
sect. 3.2.2 of Romeo (1994). 

Fig. 2 shows the contour maps of Q^^ and Qws (l^ft 
panel), and the error map of Qws (right panel). Remember 
that Q denotes the stability threshold, i.e. the value of Q* 

^ Hereafter we will use Smed and Qmedi i-e. the median values of 
s and q computed from the galaxy data of Leroy et al. (2008), 
for estimating the typical accuracy of the Wang-Silk approxima- 
tion and of our approximation. This is meant to be a complement 
to the detailed error maps shown and discussed throughout the 
paper. We do not 'hint that the stability properties can be char- 
acterized by a median value of an effective Q-parameter'. 
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Figure 2. Accuracy of the Wang-Silk approximation. The curves shown are the contour lines of Qws ^^"^ Qbr (left), and the contour lines 
of the relative error (Qws " Qbr)/Qbr (right). Q and Q denote the effective Q parameter and the stability threshold, the subscripts 
WS and BR refer to Wang & Silk (1994) and Bertin & Romeo (1988). In addition, Q* and Qg are the stellar and gaseous Toomre 
parameters, cr* and erg are the radial velocity dispersions of the two components. Note how far the Wang-Silk approximation is from the 
correct stability threshold, except for erg ~ a* (see left panel). 



above which the two-component disc is locally stable, while 
Q = Q*/ Q denotes the effective Q parameter. Both stabil- 
ity thresholds are above unity and converge to 1 + as 
s ^ 1. The first property means that a disc of stars and gas 
can be gravitationally unstable even when both components 
are separately stable, as is well known (e.g., Lin & Shu 1966). 
The second property simply means that stars and gas act as 
a single component when they have the same radial velocity 
dispersion, so that Q = KO'/7rG(E*+Sg). Apart from satisfy- 
ing those properties, the two stability thresholds arc clearly 
different. Look in particular at the contour levels 1.1-1.4 of 
Qgp^. Their slope changes abruptly across the line q — 1, 
revealing the existence of two stability phases. As discussed 
in the previous paragraph, this is an important characteris- 
tic of two-component discs, which QyfQ fails to reproduce. 
From a quantitative point of view, the error that affects 
Qws is significant but below 50% (see now the right panel 
of Fig. 2). Using the median values of s and q computed 
from the galaxy data of Leroy et al. (2008), Smed — 0.27 and 
Quied — 1-5, one finds that the typical error is about 20%. 
Remember, however, that a significant fraction of the data 
populate the two-phase region, where the error can be more 
than 40%. Note also that the error is negative, which means 
that the Wang-Silk approximation underestimates the effec- 
tive Q parameter systematically. 



2.2 Our approximation 

Let us now illustrate how to find a better approximation 
for Q. The first ingredient is to determine the asymptotic 
behaviour of Q as s — >^ and s — ^ 1. These are in fact the 
natural bounds of s. A rigorous analysis was performed by 
Romeo (1985). His results can be summarized as follows: 

(i) For s <sc 1 and q < 1, i.e. in gas-dominated stability 
regimes, Q ^ q^^ + 2s. 

(ii) For s <C 1 and g > 1, i.e. in star-dominated stability 
regimes, Q ^ 1 + 2s q^^ . 

(iii) For s « 1, i.e. in the limiting case of a one-component 
disc, Q f» 1 -h . 

Note that Q behaves asymptotically as a weighted sum of 
two terms: 1 and . Note also that the weight factors 
change symmetrically as we move from case (i) to case (iii): 

(2s, 1) — > (1,2s) — ^ (1, 1). Such symmetry suggests that we 
should search for an approximation of the form 



Q = 



W{s) + q-'^ 
l + W{s)q-^ 



if g < 1 , 
else : 



(4) 



where W{s) 2s as s — >■ 0, and W{s) « 1 as s — > 1. A 
further constraint on W{s) follows from the fact that the 
original system of fluid and Poisson equations remains unal- 
tered if we interchange the stellar and gaseous components. 
Eq. (4) must then be invariant under the transformation 
s H-^ s~^, q q~^ and hence Q i-> gQ, where qQ is the value 
of Qg above which the two-component disc is locally stable. 



4 A. B. Romeo and J. Wiegert 




Invariance requires that W{s~^) = W{s). A simple function 
that satisfies this requirement and matches the asymptotic 
behaviour above is 
2s 



W{s) = 



1 + S2 



(5) 



Since Q = Q* / Q, Eqs (4) and (5) load us to the following 
approximation for the effective Q parameter: 



1 

Q 



= < 



w_ 



1 

Qg 



— + — if > Qg 



— + — if Qg > Q* 



W ■ 



1 

Q^ 

2cr*crg 
0-2 + 0-1 



W 

Qg 



(6) 



(7) 



Our approximation is almost as simple as the Wang-Silk ap- 
proximation [Eq. (1)], but differs from that in one important 
respect: it gives less weight to the component with larger Q. 
The weight factor W depends symmetrically on the radial 
velocity dispersions of the two components, and is generally 
small. 

To evaluate the accuracy of our approximation, wc com- 
pare Q(s,q) with Qbr{s,C[) and compute the relative error 
(Q — Qbr)/Qbr as a function of s and q, as wc did for the 
Wang-Silk approximation. Fig. 3 shows that Q works well 
in the whole parameter space (see left panel). Note in par- 
ticular how successfully our approximation reproduces the 
gaseous and stellar stability phases for Q<1.2. Fig. 3 also 
shows that Q overestimates the effective stability parameter, 
but the error is well below 10% even inside the two-phase 



region (see right panel) . The error can be reduced further by 
fine-tuning the weight factor, but the approximation will no 
longer be consistent with the asymptotic behaviour of the 
stability threshold. 



2.3 How to apply our approximation to 
realistically thick discs 

As pointed out in Sect. 1, the stabilizing effect of disc thick- 
ness is usually neglected but significant. In this section, we 
show that our approximation can easily be modified so as to 
take this effect into account. 

Romeo (1992) investigated the gravitational instability 
of galactic discs taking rigorously into account two factors: 
(i) their vertical structure at equilibrium; (ii) the coupling 
between scale height, h, and vertical velocity dispersion, cr^, 
in the stellar and gaseous layers. He calculated the effective 
Q parameter both as a function of hi, and hg , and as a func- 
tion of C7z* and cr^g. He also discussed the advantages of using 
az* and cr^g as input quantities. This effective Q parameter 
has been studied further by Wiegert (2010). Hereafter we 
will denote it with Qwr,- 

Let us now illustrate how to find a simple and accurate 
approximation to Qwr. In the infinitesimally thin case [see 
Eq. (6)], the local stability level of the disc is dominated 
by the component with smaller Q. The contribution of the 
other component is weakened by the W factor, which is gen- 
erally small. This suggests that the effect of thickness can 
be estimated reasonably well by considering each component 
separately. Romeo (1994) analysed this case in detail. The 
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Figure 4. Accuracy of our approximation (left) vs. accuracy of tlie Wang-Silk approximation (right) for realistically thick discs. The 
curves shown are the contour lines of the relative errors (Q — Qwr)/Qwr (see left panel) and {Qws ~ Qwr)/Qwr (see right panel) 
for (iTz/ffl)* = 0.5 and (czl<^ii)% = 1- Here Q is our effective Q parameter, QwR is the effective Q parameter of Romeo (1992) and 
Wiegert (2010), Qws is the effective Q parameter of Wang & Silk (1994), and a^/aji is the ratio of vertical to radial velocity dispersion. 

The rest of the notation is the same as in Figs 1-3. Also shown is the corresponding two-phase region (dashed lines). The boundaries of 
this region and the transition line intersect at {s,q) ~ (0.16,0.74), where the stability threshold is QwR — 1-3. 



effect of thickness is to increase the stability parameter of 
each component by a factor T, which depends on the ratio 
of vertical to radial velocity dispersion: 



T w 0.8 -h 0.7 



(8) 



Eq. (8) can be inferred from fig. 3 (top) of Romeo (1994) 
and applies for 0.5 < a"z/o"fl ^ 1, which is the usual range of 
velocity anisotropy. To approximate Qwr, use then Eq. (6) 
with Qt, and Qg replaced by T^Q* and TgQg-. 



1 
Q 



( W 1 

+ if T,Q, > TgQg , 



1 w 

-I- -^^TTT if TgQg > ; 



(9) 



k T*(5* TgQg 



where Q is our effective Q parameter for realistically thick 

discs, W is given by Eq. (7), and Tg are given by Eq. 
(8). Eq. (9) tells us that the local stability level of the disc 
is now dominated by the component with smaller TQ. The 
contribution of the other component is still suppressed by 
the W factor. 

The left panel of Fig. 4 shows the error map of Q for a 
galactic disc with (az/cFR)-/, = 0.5 and (az/crB)g = 1, and the 
corresponding two-phase region (Wiegert 2010). Note that 
the error is below 15% even inside this region, which con- 
firms the high accuracy of our approximation in this more 
realistic context. What about the accuracy of the Wang-Silk 
approximation? The right panel of Fig. 4 shows that the rel- 
ative error (Qws— Qwr)/Qwr is much larger than ours, and 
can be well above 50% inside the two-phase region. 



2.4 Application to nearby star-forming spirals 

In this section, we show how to use our effective Q param- 
eter for measuring the stability level of galactic discs, and 
why such a diagnostic is more predictive than the classical 
Toomre parameter. 

We consider the same sample of spiral galaxies as in 
Sect. 2.1, and refer to Leroy et al. (2008) for a detailed 
description of the data and their translation into physical 
quantities. For each galaxy, we compute the radial profile 
of our effective Q parameter, Q, using Eq. (9). We adopt 
(o'z/o'r)* = 0.6, as was assumed by Leroy et al. (2008), 
and {(7z/o'B)g ~ 1, as is natural for a collisional component. 
We also compute the radial profile of the gaseous Toomre 
parameter, Qg, which is the traditional diagnostic used for 
predicting star formation thresholds in galaxies (e.g.. Quirk 
1972; Kennicutt 1989; Martin & Kennicutt 2001; Schaye 
2008; Elmegreen 2011). 

Fig. 5 shows Qg{R) and Q{R) for the whole galaxy sam- 
ple. Note that Qg spans a much wider range of values than 
Q at any given R. This is true even at distances as large 
as the optical radius, R2r,, where Qg is supposed to be a 
reliable diagnostic. A similar fact was noted by Leroy et 
al. (2008), using an effective Q parameter that neglects the 
stabilizing effect of disc thickness (Jog 1996; Rafikov 2001). 
Why are Qg and Q so weakly correlated across the entire 
optical disc? Eq. (9) helps us to clarify this point. It tells 
us that the value of Q is dominated by the gaseous com- 
ponent if TgQg < TtQi,, and by the stellar component if 
T*Q* < TgQg. In the right panel of Fig. 5, we have colour- 
coded the data so as to show whether TgQg < T*Q« or vice 
versa. It turns out that in 92% of the cases the value of Q 
is dominated by the stellar component. Gas dominates the 
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Figure 5. The stability level of nearby star-forming spirals, as measured by two diagnostics: the gaseous Toomre parameter, Qg, and 
our effective Q parameter, Q [see Eq. (9)]. The galaxy data are from Leroy et al. (2008), R is the galactocentric distance, and R25 is 
the optical radius. In the right panel, the data are colour-coded so as to show whether the stability level is gas- or star- dominated, as 
predicted by Eq. (9). The two data points that lie well below the critical stability level tell us that the nuclear region of NGC 6946 is 
subject to strong gas-dominated instabilities. This is consistent with the fact that NGC 6946 hosts a nuclear starburst (e.g., Engelbracht 
et al. 1996). 



stability level only in 8% of the cases. This is why Qg and 
Q are so weakly correlated. This result illustrates (i) how 
important it is to consider both gas and stars when measur- 
ing the stability level of galactic discs, and (ii) the strong 
advantage of using our effective Q parameter as a stability 
diagnostic. 



3 CONCLUSIONS 

• The approximation of Wang & Silk (1994) [Eq. (1)] un- 
derestimates the effective Q parameter. The error is typically 
20%, but can be as large as 40% or more if Ug <0.2(j* and 
Qg ~ Q*. In this case, the gaseous and stellar components 
should contribute separately to the gravitational instability 
of the disc (Berlin & Romeo 1988). But such dynamical de- 
coupling is difficult to approximate because it involves two 
stability regimes, one dominated by the gas and the other 
dominated by the stars, and because there is a sharp tran- 
sition between the two 'phases'. So it is not strange that 
the Wang-Silk approximation becomes less accurate when 
(Tg < 0.2(T* and Qg ~ Q*. 

• Our approximation [Eq. (6)] overestimates the effective 
Q parameter, but the error is less than 9% and typically as 
small as 4%. The accuracy and simplicity of our approxima- 
tion result from a rigorous analysis, which takes into account 
the stability characteristics of the disc as well as the sym- 
metries of the problem. 

• We provide a simple recipe for applying our approxi- 
mation to realistically thick discs [see Eq. (9)]. The ratio of 
vertical to radial velocity dispersion is usually 0.5 for the 
stars and 1 for the gas. In this case, our approximation is in 
error by less than 15%, whereas the Wang-Silk approxima- 
tion can be in error by more than 50%. Note also that the 
effective Q parameter is 20-50% larger than in the infinitesi- 
mally thin case. Thus the effect of thickness is important and 



should be taken into account when analysing the stability 
of galactic discs. 
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